function [wiener_result,H]=my_wiener(noisy,clean,noise,para)%加噪信号；干净信号；噪声

    alpha=para.alpha;
    beta=para.beta;
    
    s_noisy=stft(noisy);
    s_noise=stft(noise);
    s_clean=stft(clean);
    
    pxx=mean((abs(s_clean)).^2);%模值平方的期望
    pnn=mean((abs(s_noise)).^2);
    
    H=(pxx./(pxx+alpha*pnn)).^beta;
    
    s_enhec=s_noisy.*H;
    wiener_result=istft(s_enhec);
end